Method for determining the humidity and density of a dielectric material 

The invention relates to a method for determining the humidity and/or density 
of a dielectric material in a resonator filled with the material, comprising a 
sender and a receiver. 

The dielectric properties of a material, described by the complex relative 
dielectric constant fir = e[ - icr, can be affected by humidity and density in case 
of porous materials. Humidity and density thus modify the scalar parameters 
resonant frequency fr and resonator quality Q of a material-filled resonator in 

contrast to those of the empty, air-filled resonator (/^o Qo) in such a way 
that 

" and Qm Qo , wherein 1 < fir < eL and 0 < < eJL, wherein 
and ^rLx are the maximum values resulting for the respective material 
from the assigned humidity range and density range. 

In prior art, different methods are known for measuring the humidity or density 
of granular materials, in which the resonance behavior of a matter-filled 
resonator is used. 

For example, from US 566 60 61 a method is known for measuring the 
humidity in granular materials by means of microwaves. There, one 
electronically reacts to a threshold while bidirectionally sweeping a frequency 
range. From the chronological properties of pulses created in such a way the 
components of the dielectric constant are deduced. 

All methods are relatively slow and have either a dependency of the density of 
the material or have relatively large errors in the determination of the humidity. 
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The invention is underlied by the problem to specify a method of the type 
initially mentioned, by which a fast, accurate and density-independent 
determination of the humidity is possible. 

According to the invention, the problem is solved by a method comprising the 
attributes given in claim 1. 

Advantageous embodiments are given in the dependent claims. 

For the purpose of the invention, any arbitrary quantity describing the 
resonance width can be considered as a bandwith in the following, wherein 
appropriate adjustments have to be provided in the respective definitions, in 
particular of the threshold values, and equations. 

By digitally recording the resonance curve a fast acquisition of the measuring 
values is possible. Thus a close chronological-spatial assignment of the 
humidity determined from these values to material dynamically guided through 
the resonator is possible. For this purpose, the invention provides for that the 
sender emits a signal; a resonance curve of the filled resonator is sweeped, 
wherein respective related signal strength values of the receiver signal are 
measured at different frequencies; the resonant frequency and the bandwidth 
are determined for the filled resonator from the measured points; and the 
humidity and/or density of the material is calculated by solving a second 
system of equations comprising the resonant frequencies and the bandwidths of 
the empty and of the filled resonator and known calibration coefficients of the 
resonator. This method enables a density-independent determination of the 
humidity of the material. Besides, the density of the material can be acquired 
with little effort. 



A preferred embodiment provides for that, from the points for determining the 
bandwith of the filled resonator, either the quantities resonant frequency, 
resonator quality and resonance maximum are determined and the bandwidth is 
calculated therefrom, or cut-off frequencies are determined and the resonant 
5 frequency and the bandwidth are calculated therefrom. 

In another embodiment a lower threshold value is calculated and a second 
sweeping pass with smaller step sizes is performed in that range in which the 
signal strength values are higher than the threshold value. By a two-pass 
10 procedure the accuracy can be significantly increased and, nevertheless, the 

required time can be kept short, in particular if, during the second sweeping 
pass, a decreased step size is used in a range of the resonance peak only. 

Advantageously, sweeping the resonance curve is performed in equally spaced 
15 steps. The simplest and fastest form of sweeping consists of equally spaced 

steps. Variable step sizes in the second sweeping pass can shorten its duration. 

Preferably, the sender is operated using a constant strength. The measuring 
values at the receiver can be used without adjusting or scaling them in case of a 
20 constant signal strength. 

In a possible embodiment the cut-off frequencies of the resonator are 
determined by determining the point having the highest receiver signal strength 
value, and, starting from this point, calculating a threshold value; determining 
25 two respective proximate points for positive and negative slope sections, the 

signal strength values of these points lying below and above the threshold 
value; calculating first and second cut-off frequencies therefrom by 
respectively interpolating between the proximate points. Determining the cut- 
off frequencies by interpolating between point pairs surrounding a threshold 
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value as an initial parameter for determining the humidity and/or density is a 
fast and simple method. 

An advantageous embodiment thereby provides for that the threshold value 
corresponds to an attenuation of 3 dB in relation to the highest signal value. If 
the threshold value is chosen corresponding to an attenuation of 3 dB, starting 
from the maximum signal strength value, the equations to be solved obtain a 
very simple form. 

In an alternative embodiment the quantities resonant frequency, resonator 
quality and resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting three points and solving a first system of equations 
for these quantities, the system consisting of three equations of an analytic 
resonance curve that are valid for the three points. The resonant frequency, the 
resonator quality and the resonance ampitude can be directly determined as the 
initial parameters for determining the humidity and/or density even faster and 
with less error than by interpolation by solving a first system of equations 
consisting of three equations and thus being completely determined, wherein 
for each point out of a group of three one of the three equations is valid, the 
three points being selected from the present set. In particular, noisy resonance 
curves can be analyzed more exactly this way. 

In another alternative embodiment the quantities resonant frequency, resonator 
quality, resonance maximum of the resonator are determined by arbitrarily 
and/or randomly selecting a set of points whose number is an integer multiple 
of three and at least six, and splitting up the point set into three equally sized 
groups; for each combination of three points, wherein each point comes from a 
different group, solving a first system of equations for these quantities, the 
system consisting of three equations of the analytic resonance curve valid for 



these three points; and creating the average for each quantity from the values 
calculated at the combinations. Even more exact values than with three points 
can be obtained by creating several triple-groups from the present point set and 
averaging the initial parameters obtained thereby over all groups. 

In an advantageous embodiment, as a condition for arbitrarily and/or randomly 
selecting the points, the signal value of a point to be selected is higher than the 
highest signal value attenuated by 3 dB. In both procedures using three points 
or a multiple thereof, preferably only those points from the resonance curve are 
selected whose signal strength is higher than the maximum value of all 
measurement values, attenuated by 3 dB, as this way only significant values are 
used. 

The system of equations to be solved is preferably chosen in such a way that it 
describes, in a good approximation, the correlation of humidity and density 
with a variation of resonant frequency and resonator quality or, respectively, 
with a variation of resonant frequency and bandwidth. 

Advantageously the second system of equations is non-linear. 

In a preferred embodiment the sweeping by means of the sender is performed 
up to the microwave area. 

Advantageously, voltage values or current values of the receiver are used for 
measuring the receiver signal Preferably, the electrical voltage rising at the 
receiver serves for acquiring the resonance curve, because it is measurable 
easily and without back coupling. However, the current in a receiver circuit can 
be measured, too. 



In the following, the denomination Ui is used for the measured quantities. 
However, in doing so, it shall not be implied to use the voltage only. 

The resonance curve to be reconstructed from the discrete measuring points is 
an approximation as-good-as-possible to the real resonance curve only. It is 
characterized by the course of the signal amplitude U at the receiver against the 

supplied frequency wherein Ur is the resonance maximum and Q is the 
quality of the resonance: 

The cut-off frequencies fa and fb are those frequencies at which the signal 
strength exceeds and falls below a defined value, respectively. For this purpose, 
preferably the value corresponding to an attenuation of the maximum by 3 dB 

is chosen: ^ ^ JI , Between the cut-off frequencies and the curve 
parameters resonant frequency fr and resonator quality Q there are the 

correlations: VfW7 , Q = The distance between the cut-off 

frequencies is defined as the bandwidth = fb - fa = -q 



In the following, the invention is described in further detail using examples of 
20 embodiments. 

For this purpose, 

Figure 1 shows a sweeped resonance curve, 
Figure 2 shows a resonance curve sweeped in a two-pass procedure, 
25 Figure 3 shows a first way for determining the resonance parameters, 

Figure 4 shows a second way for determining the resonance parameters. 
Figure 5 shows a third way for determining the resonance parameters 
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and 

Figure 6 shows a two-dimensional representation of calibration curves. 



In Figure 1, two resonance curves are depicted above each other like they have 
5 been recorded using the method according to the invention; the right resonance 

curve at an empty, the left resonance curve at a material-filled resonator. 

By sweeping (wobbling) the resonator between the start frequency fstarti and 
the stop frequency ^stopi the resonance curve is acquired in discrete steps. The 
10 start frequency is calculated from the maximally shifted fi-equency 

f J i_ ff 

■^"^ V^^max , the maximally changed quality Qm " Qo ^rmax and the normalized 

voltage ratio t^^" ^ : ^stani - 2Q„ V "♦"V^zq^J • +ir™ jYie 
stop frequency, using f^r = fro, results as: 

f - f»o / , // fiO \? l-a^ . ro 

15 

The sweeping speed depends, among others, on the number "i of sweeping 
points and can be increased by a two-pass procedure as depicted in Figure 2. 
For this purpose, first a sweeping is performed using a smaller number of 

A-P fstopl~fstaftl 

sweeping points and, hence, a larger frequency step size ~ ni-i and, in 
20 a second iteration, after previously determining the start and stop frequencies 

fstara =/| U«Ua Ar>0, fstopl =/| U^Ua A/<0 or fstara = f\ U^Ua A U<Vr, 

fstop2=f\iMUaAU>Ur^ another sweeping is subsequently performed using a 

A f fstop2~fstart2 

smaller frequency step size nz-i between these frequencies . 

25 From the sweeped, measured resonance curve the resonator parameters 

fr,QandUr can be determined, for example, according to the fafb procedure, 
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the three-points procedure or the 3k-points procedure. Other procedures are 
possible, too. 

Figure 3 shows the fafb procedure. It is based on directly determining a first 
and a second cut-off frequency fa, fb from the measured resonance curve. For 
this purpose, first the point (also called sweeping point) 
(f = fmax/U = Umax)having the highest voltage U^Umax is determined for 

calculating a 3 dB threshold line ^''^ ~ 72" . Subsequently, the respective two 
points in the immediate proximity of the 3 dB threshold line are determined. By 
linearly interpolating between the sweeping points al and a2 as well as bl and 
b2 the first and the second 3 dB cut-off frequency fa, fb are obtained: 

fa " fal + "U;^-U7"(fa2 " fal)^ fb = fbl + Ub2-Ubi (fb2 " fbl). 

From the equations ^ T^Wb , Q = fj-f, ^ fr and Q can be calculated. The 
resonance voltage Ur then results as: 

Calculating the cut-off frequencies is more or less error-prone because of the 
discretization of the resonance curve - on one hand, by the determination of the 
3 dB threshold line (Ua = < Ur/V2 ) 

derived from Umax and, on the 
other hand, by the linear interpolation between the sweeping points al, a2, bl, 
b2. A higher sweeping rate indeed reduces the errors, but also increases the 
sweeping time and thus reduces the sweeping speed. At best, a non-linear 
interpolation is possible in this procedure in order to improve accuracy. 
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In the three-points procedure shown in Figure 4, three sweeping points 
U,,U2,U3 are selected arbitrarily or by a random generator preferably above 
the line Us = s • Umax with ^^1^. 



By solving the first system of equations Gl: 



the resonator parameters/'-' Q Ur obtained: 




I u?y) , u^(u|-u|)^u|(i| 




The advantage of the three-points procedure in relation to the fafb procedure 
lies in that no errors occur because of the discretization of the resonance curve, 
and that only few points are needed for sweeping the resonance curve and thus 
a high sweeping speed is achieved. 

If the resonance curve is noisy, error-prone resonator parameters are 
determined by the three-points procedure, too. hi order to reduce the impact of 
the noise more than three sweeping points are selected, namely a 
number^'k, k = 2,3,4,.,.^ and are split up into three groups Ml, M2, M3, 
having k points per group above the line Us = s • Umax . Figure 5 shows this 3k- 
points procedure. For all k^ possibilities of combining one respective point 
from each of the three groups, the resonator parameters are determined from 
the first system of equations Gl according to the three-points procedure 

described above, and the k^ partial resuhs are averaged subsequently: 
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In order to determine humidity and density from determined resonator 
parameters the used arrangement is calibrated using material of known 
5 humidity and density. The calibration is advantageously performed out of the 

regular operation. 

In the calibration, the electrical resonator values resonant frequency fr, 
resonator quality Q and bandwidth BW(BW = fr/Q) are respectively assigned 
to the material quantities humidity ¥ and density Q, This assignment between 
the material values and the electrical values is non-linear and can be described 
for a defined range of humidity and density with sufficient accuracy, for 
example, by the following second system of equations G2: 
^fr = ^friQ^ + af,, ^ bf,2^2 \i/ + hf^,Q\i/ + Cf,,Q^ + Cf,, g 
ABW = aewzp^ + ^bw, ev^^ + bewz^^ y/ + bew, + cbw2^^ + caw, Q 

^fr and ABW are the differences of resonant frequency and bandwidth, 
respectively, between the empty and the material-filled resonator: 

Afr = fro - fr., ABW = BWn, - BWq 

Now, it is an object of the calibration to determine the twelve calibration 
coefficients ^f^y bf„, bf,2, Cf,,, Cf,2» ^bw,, aswz, baw,, bawj, cawi, cawz 
for the respective material from a sufficient large number of calibration values 
(Afr and ABW ^^ith related material values W und Q). 

For this purpose, the measuring values of the resonant frequency difference Afr 
and the bandwidth difference ABW are assigned to the humidity and density 
values that have been determined by appropriate reference methods. These 
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calibration values characterizing the resonant frequency and bandwidth 
depending on humidity and density are the basis for determining the calibration 
coefficients. The regression of the calibration values of equal humidity is 
performed in the form of Afr = df,2e^ + df,i^ and Afr = dBW2^^ + dBw,^, 
wherein the regression curves have to run through the origin, as and ABW 
are zero for the empty resonator (air having 6^^). The regression provides 
calibration points of equal density depending on the humidity. Using these, 
another regression of the form ^fr = ar, + bf, + Cf, and 
ABW= aewy^^+bewV^ + CBw is performed. The values determined this way 
for the regression coefficients ^fr» bf,, Cf, and ^bw, bew, Cbw are plotted 
against the density and, therefrom, the calibration coefficients are determined 
by a quadratic regression. 

Figur 6 shows the results of the calibratrion in the form of a set of curves which 
also serves for analysis. 

From the measured resonant frequency difference values and bandwidth 
difference values W humidity content ¥ and density 6 are calculated for the 
respective material by solving the above second system of equations G2. 

In doing so, two real and two imaginary roots do result. From the course of the 
calibration curves in the Afr.ABW diagram it can be determined if there is only 
one real solution in the humidity and density range of interest. In the Afr_ABW 
diagram the bandwidth difference is represented in dependency of the resonant 
frequency difference for curves of equal density and humidity. If the course of 
these curves, characterized by the points A, B, C und D in Figure 6, is 
continuous and imique in the himiidity and density range of interest then only 
one real solution exists in this range. 
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For solving the second system of equations G2 an iterative procedure is 
appropriate. For this purpose, the second system of equations G2 is solved for 

y - 2afr2C^2afri(? ^ l^2af,2(?+2afri(?>/ ^U2Q'^-^^U\Q 

_ bBW2(>+bBWi //^ bBW2(?+bBWi CBW2g^+CBWi-ABWo 

~ 2aBW2(?+2aBWie ^ V.2aBW2^2aBW|g J "~ aBW2(?^+aBW|(? 

From the intersection of both equations in a humidity-density diagram the 
sought values for ¥ and Q are obtained. 
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List of reference numerals 



fr 

fro 

fr. 

f Start 1 > fstopl 
Umax 

a, S 

fstart2j fstop2 
fmax 

(fal/Ual), (fa2/Ua2) 
(fbl/Ubl), (fb2/Ub2) 
(fl/U,)... (f3/U3) 

M1,M2, M3 
(fu/Uii)...(fi4/U,4) 

(f2l/U2l).,.(f24/U24) 
(f3,/U3l) ... (f34/U34) 

A,B,C,D 



Resonant frequency in the general case 

Resonant frequency of the empty resonator 

Resonant frequency of the filled resonator 

Start and stop frequencies for first sweeping pass 

Highest signal strength value measured 
Threshold value factors 

Start and stop frequencies for second sweeping pass 

Frequency at which the highest signal strength 
value is present 

Proximate points of the first cut-off frequency 

Proximate points of the second cut-off frequency 

Three selected points 
point groups 

Elements of point group Ml 

Elements of point group M2 

Elements of point group M3 

Limits of the humidty and density range of interest 
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